Representation and quantification of module activity from omics data with rROMA.

The efficiency of analyzing high-throughput data in systems biology has been demonstrated in numerous studies, where molecular data, such as transcriptomics and proteomics, offers great opportunities for understanding the complexity of biological processes. One important aspect of data analysis in systems biology is the shift from a reductionist approach that focuses on individual components to a more integrative perspective that considers the system as a whole, where the emphasis shifted from differential expression of individual genes to determining the activity of gene sets. Here, we present the rROMA software package for fast and accurate computation of the activity of gene sets with coordinated expression. The rROMA package incorporates significant improvements in the calculation algorithm, along with the implementation of several functions for statistical analysis and visualizing results. These additions greatly expand the package's capabilities and offer valuable tools for data analysis and interpretation. It is an open-source package available on github at: www.github.com/sysbio-curie/rROMA . Based on publicly available transcriptomic datasets, we applied rROMA to cystic fibrosis, highlighting biological mechanisms potentially involved in the establishment and progression of the disease and the associated genes. Results indicate that rROMA can detect disease-related active signaling pathways using transcriptomic and proteomic data. The results notably identified a significant mechanism relevant to cystic fibrosis, raised awareness of a possible bias related to cell culture, and uncovered an intriguing gene that warrants further investigation.

A Network of 17 Microtubule-Related Genes Highlights Functional Deregulations in Breast Cancer.

A wide panel of microtubule-associated proteins and kinases is involved in coordinated regulation of the microtubule cytoskeleton and may thus represent valuable molecular markers contributing to major cellular pathways deregulated in cancer. We previously identified a panel of 17 microtubule-related (MT-Rel) genes that are differentially expressed in breast tumors showing resistance to taxane-based chemotherapy. In the present study, we evaluated the expression, prognostic value and functional impact of these genes in breast cancer. We show that 14 MT-Rel genes (KIF4A, ASPM, KIF20A, KIF14, TPX2, KIF18B, KIFC1, AURKB, KIF2C, GTSE1, KIF15, KIF11, RACGAP1, STMN1) are up-regulated in breast tumors compared with adjacent normal tissue. Six of them (KIF4A, ASPM, KIF20A, KIF14, TPX2, KIF18B) are overexpressed by more than 10-fold in tumor samples and four of them (KIF11, AURKB, TPX2 and KIFC1) are essential for cell survival. Overexpression of all 14 genes, and underexpression of 3 other MT-Rel genes (MAST4, MAPT and MTUS1) are associated with poor breast cancer patient survival. A Systems Biology approach highlighted three major functional networks connecting the 17 MT-Rel genes and their partners, which are centered on spindle assembly, chromosome segregation and cytokinesis. Our studies identified mitotic Aurora kinases and their substrates as major targets for therapeutic approaches against breast cancer.

Exploring isofunctional molecules: Design of a benchmark and evaluation of prediction performance.

Identification of novel chemotypes with biological activity similar to a known active molecule is an important challenge in drug discovery called 'scaffold hopping'. Small-, medium-, and large-step scaffold hopping efforts may lead to increasing degrees of chemical structure novelty with respect to the parent compound. In the present paper, we focus on the problem of large-step scaffold hopping. We assembled a high quality and well characterized dataset of scaffold hopping examples comprising pairs of active molecules and including a variety of protein targets. This dataset was used to build a benchmark corresponding to the setting of real-life applications: one active molecule is known, and the second active is searched among a set of decoys chosen in a way to avoid statistical bias. This allowed us to evaluate the performance of computational methods for solving large-step scaffold hopping problems. In particular, we assessed how difficult these problems are, particularly for classical 2D and 3D ligand-based methods. We also showed that a machine-learning chemogenomic algorithm outperforms classical methods and we provided some useful hints for future improvements.

Differential CFTR-Interactome Proximity Labeling Procedures Identify Enrichment in Multiple SLC Transporters.

Proteins interacting with CFTR and its mutants have been intensively studied using different experimental approaches. These studies provided information on the cellular processes leading to proper protein folding, routing to the plasma membrane, recycling, activation and degradation. Recently, new approaches have been developed based on the proximity labeling of protein partners or proteins in close vicinity and their subsequent identification by mass spectrometry. In this study, we evaluated TurboID- and APEX2-based proximity labeling of WT CFTR and compared the obtained data to those reported in databases. The CFTR-WT interactome was then compared to that of two CFTR (G551D and W1282X) mutants and the structurally unrelated potassium channel KCNK3. The two proximity labeling approaches identified both known and additional CFTR protein partners, including multiple SLC transporters. Proximity labeling approaches provided a more comprehensive picture of the CFTR interactome and improved our knowledge of the CFTR environment.

Profiling the response to lumacaftor-ivacaftor in children with cystic between fibrosis and new insight from a French-Italian real-life cohort.

INTRODUCTION: Clinical trials for CFTR modulators consider mean changes of clinical status at the cohort level, and thus fail to assess the heterogeneity of the response. We aimed to study the different response profiles to lumacaftor-ivacaftor according to age in children with cystic fibrosis (CF). METHODS: A mathematical framework, including principal component analysis, data clustering, and data completion, was applied to a multicenter cohort of 112 children aged 6-18 years, treated with lumacaftor-ivacaftor. Studied parameters at baseline and 6 months included body mass index (BMI), number of days of antibiotics (ATB), Sweat test (ST), forced expiratory volume in 1 s expressed in percentage predicted (ppFEV1 ), forced vital capacity (ppFVC), and forced expiratory flow at 25%-75% of FVC (ppFEF25-75 ). RESULTS: Change in ppFEV1 was the most significant parameter in characterizing response heterogeneity among the 12-18-year-old patients. Patients with minimal changes in ppFEV1 were further separated by change in BMI and ATB course. In the 6-12-year-old children both BMI and ppFEV1 evolution were the most relevant. ST change was not associated with a clinical response. CONCLUSIONS: Change in ppFEV1 , BMI, and ATB course are the most relevant outcomes to discriminate clinical response profiles in children treated with lumacaftor-ivacaftor. Prepubertal and pubertal children display different response profiles.

Drug Target Identification with Machine Learning: How to Choose Negative Examples.

Identification of the protein targets of hit molecules is essential in the drug discovery process. Target prediction with machine learning algorithms can help accelerate this search, limiting the number of required experiments. However, Drug-Target Interactions databases used for training present high statistical bias, leading to a high number of false positives, thus increasing time and cost of experimental validation campaigns. To minimize the number of false positives among predicted targets, we propose a new scheme for choosing negative examples, so that each protein and each drug appears an equal number of times in positive and negative examples. We artificially reproduce the process of target identification for three specific drugs, and more globally for 200 approved drugs. For the detailed three drug examples, and for the larger set of 200 drugs, training with the proposed scheme for the choice of negative examples improved target prediction results: the average number of false positives among the top ranked predicted targets decreased, and overall, the rank of the true targets was improved.Our method corrects databases' statistical bias and reduces the number of false positive predictions, and therefore the number of useless experiments potentially undertaken.

Urinary Exosomes of Patients with Cystic Fibrosis Unravel CFTR-Related Renal Disease.

Background: The prevalence of chronic kidney disease is increased in patients with cystic fibrosis (CF). The study of urinary exosomal proteins might provide insight into the pathophysiology of CF kidney disease. Methods: Urine samples were collected from 19 CF patients (among those 7 were treated by cystic fibrosis transmembrane conductance regulator (CFTR) modulators), and 8 healthy subjects. Urine exosomal protein content was determined by high resolution mass spectrometry. Results: A heatmap of the differentially expressed proteins in urinary exosomes showed a clear separation between control and CF patients. Seventeen proteins were upregulated in CF patients (including epidermal growth factor receptor (EGFR); proteasome subunit beta type-6, transglutaminases, caspase 14) and 118 were downregulated (including glutathione S-transferases, superoxide dismutase, klotho, endosomal sorting complex required for transport, and matrisome proteins). Gene set enrichment analysis revealed 20 gene sets upregulated and 74 downregulated. Treatment with CFTR modulators yielded no significant modification of the proteomic content. These results highlight that CF kidney cells adapt to the CFTR defect by upregulating proteasome activity and that autophagy and endosomal targeting are impaired. Increased expression of EGFR and decreased expression of klotho and matrisome might play a central role in this CF kidney signature by inducing oxidation, inflammation, accelerated senescence, and abnormal tissue repair. Conclusions: Our study unravels novel insights into consequences of CFTR dysfunction in the urinary tract, some of which may have clinical and therapeutic implications.

Evaluation of deep and shallow learning methods in chemogenomics for the prediction of drugs specificity.

Chemogenomics, also called proteochemometrics, covers a range of computational methods that can be used to predict protein-ligand interactions at large scales in the protein and chemical spaces. They differ from more classical ligand-based methods (also called QSAR) that predict ligands for a given protein receptor. In the context of drug discovery process, chemogenomics allows to tackle the question of predicting off-target proteins for drug candidates, one of the main causes of undesirable side-effects and failure within drugs development processes. The present study compares shallow and deep machine-learning approaches for chemogenomics, and explores data augmentation techniques for deep learning algorithms in chemogenomics. Shallow machine-learning algorithms rely on expert-based chemical and protein descriptors, while recent developments in deep learning algorithms enable to learn abstract numerical representations of molecular graphs and protein sequences, in order to optimise the performance of the prediction task. We first propose a formulation of chemogenomics with deep learning, called the chemogenomic neural network (CN), as a feed-forward neural network taking as input the combination of molecule and protein representations learnt by molecular graph and protein sequence encoders. We show that, on large datasets, the deep learning CN model outperforms state-of-the-art shallow methods, and competes with deep methods with expert-based descriptors. However, on small datasets, shallow methods present better prediction performance than deep learning methods. Then, we evaluate data augmentation techniques, namely multi-view and transfer learning, to improve the prediction performance of the chemogenomic neural network. We conclude that a promising research direction is to integrate heterogeneous sources of data such as auxiliary tasks for which large datasets are available, or independently, multiple molecule and protein attribute views.

LOTUS: A single- and multitask machine learning algorithm for the prediction of cancer driver genes.

Cancer driver genes, i.e., oncogenes and tumor suppressor genes, are involved in the acquisition of important functions in tumors, providing a selective growth advantage, allowing uncontrolled proliferation and avoiding apoptosis. It is therefore important to identify these driver genes, both for the fundamental understanding of cancer and to help finding new therapeutic targets or biomarkers. Although the most frequently mutated driver genes have been identified, it is believed that many more remain to be discovered, particularly for driver genes specific to some cancer types. In this paper, we propose a new computational method called LOTUS to predict new driver genes. LOTUS is a machine-learning based approach which allows to integrate various types of data in a versatile manner, including information about gene mutations and protein-protein interactions. In addition, LOTUS can predict cancer driver genes in a pan-cancer setting as well as for specific cancer types, using a multitask learning strategy to share information across cancer types. We empirically show that LOTUS outperforms five other state-of-the-art driver gene prediction methods, both in terms of intrinsic consistency and prediction accuracy, and provide predictions of new cancer genes across many cancer types.

Efficient multi-task chemogenomics for drug specificity prediction.

Adverse drug reactions, also called side effects, range from mild to fatal clinical events and significantly affect the quality of care. Among other causes, side effects occur when drugs bind to proteins other than their intended target. As experimentally testing drug specificity against the entire proteome is out of reach, we investigate the application of chemogenomics approaches. We formulate the study of drug specificity as a problem of predicting interactions between drugs and proteins at the proteome scale. We build several benchmark datasets, and propose NN-MT, a multi-task Support Vector Machine (SVM) algorithm that is trained on a limited number of data points, in order to solve the computational issues or proteome-wide SVM for chemogenomics. We compare NN-MT to different state-of-the-art methods, and show that its prediction performances are similar or better, at an efficient calculation cost. Compared to its competitors, the proposed method is particularly efficient to predict (protein, ligand) interactions in the difficult double-orphan case, i.e. when no interactions are previously known for the protein nor for the ligand. The NN-MT algorithm appears to be a good default method providing state-of-the-art or better performances, in a wide range of prediction scenario that are considered in the present study: proteome-wide prediction, protein family prediction, test (protein, ligand) pairs dissimilar to pairs in the train set, and orphan cases.

Design, synthesis, biological evaluation and cellular imaging of imidazo[4,5-b]pyridine derivatives as potent and selective TAM inhibitors.

The TAM kinase family arises as a new effective and attractive therapeutic target for cancer therapy, autoimmune and viral diseases. A series of 2,6-disubstituted imidazo[4,5-b]pyridines were designed, synthesized and identified as highly potent TAM inhibitors. Despite remarkable structural similarities within the TAM family, compounds 28 and 25 demonstrated high activity and selectivity in vitro against AXL and MER, with IC50 value of 0.77 nM and 9 nM respectively and a 120- to 900-fold selectivity. We also observed an unexpected nuclear localization for compound 10Bb, thanks to nanoSIMS technology, which could be correlated to the absence of cytotoxicity on three different cancer cell lines being sensitive to TAM inhibition.

Kernel Multitask Regression for Toxicogenetics.

The development of high-throughput in vitro assays to study quantitatively the toxicity of chemical compounds on genetically characterized human-derived cell lines paves the way to predictive toxicogenetics, where one would be able to predict the toxicity of any particular compound on any particular individual. In this paper we present a machine learning-based approach for that purpose, kernel multitask regression (KMR), which combines chemical characterizations of molecular compounds with genetic and transcriptomic characterizations of cell lines to predict the toxicity of a given compound on a given cell line. We demonstrate the relevance of the method on the recent DREAM8 Toxicogenetics challenge, where it ranked among the best state-of-the-art models, and discuss the importance of choosing good descriptors for cell lines and chemicals.

MULTITASK FEATURE SELECTION WITH TASK DESCRIPTORS.

Machine learning applications in precision medicine are severely limited by the scarcity of data to learn from. Indeed, training data often contains many more features than samples. To alleviate the resulting statistical issues, the multitask learning framework proposes to learn different but related tasks jointly, rather than independently, by sharing information between these tasks. Within this framework, the joint regularization of model parameters results in models with few non-zero coefficients and that share similar sparsity patterns. We propose a new regularized multitask approach that incorporates task descriptors, hence modulating the amount of information shared between tasks according to their similarity. We show on simulated data that this method outperforms other multitask feature selection approaches, particularly in the case of scarce data. In addition, we demonstrate on peptide MHC-I binding data the ability of the proposed approach to make predictions for new tasks for which no training data is available.

Epigenomic alterations in breast carcinoma from primary tumor to locoregional recurrences.

INTRODUCTION: Epigenetic modifications such as aberrant DNA methylation has long been associated with tumorogenesis. Little is known, however, about how these modifications appear in cancer progression. Comparing the methylome of breast carcinomas and locoregional evolutions could shed light on this process. METHODS: The methylome profiles of 48 primary breast carcinomas (PT) and their matched axillary metastases (PT/AM pairs, 20 cases), local recurrences (PT/LR pairs, 17 cases) or contralateral breast carcinomas (PT/CL pairs, 11 cases) were analyzed. Univariate and multivariate analyzes were performed to determine differentially methylated probes (DMPs), and a similarity score was defined to compare methylation profiles. Correlation with copy-number based score was calculated and metastatic-free survival was compared between methods. RESULTS: 49 DMPs were found for the PT/AM set, but none for the others (FDR < 5%). Hierarchical clustering clustered 75% of the PT/AM, 47% of the PT/LR, and none of the PT/CL pairs together. A methylation-based score (MS) was defined as a clonality measure. The PT/AM set contained a high proportion of clonal pairs while PT/LR pairs were evenly split between high and low MS score, suggesting two groups: true recurrences (TR) and new primary tumors (NP). CL were classified as new tumors. MS score was significantly correlated with copy-number based scores. There was no significant difference between the metastatic-free survival of groups of patients based on different classifications. CONCLUSION: Epigenomic alterations are well suited to study clonality and track cancer progression. Methylation-based classification of TR and NP performed as well as clinical and copy-number based methods suggesting that these phenomenons are tightly linked.

New aminopyrimidine derivatives as inhibitors of the TAM family.

In this study, we describe the synthesis of new pyrimidine analogs of BMS-777607, a potent and selective inhibitor of Met kinase. Inhibition of Met and Axl remained high whereas inhibition of Tyro3 and Mer decreased to some extend. The preferential moderate inhibition of the non-phosphorylated form of Abl1 of some derivatives suggests that they behave as type II inhibitors. This hypothesis was confirmed by docking studies into the structure of Met (3F82) and in a Tyro3 model where key interactions with the hinge region, the DFG-out motif and the allosteric pocket explain this inhibition.

In Silico screening on the three-dimensional model of the Plasmodium vivax SUB1 protease leads to the validation of a novel anti-parasite compound.

Widespread drug resistance calls for the urgent development of new antimalarials that target novel steps in the life cycle of Plasmodium falciparum and Plasmodium vivax. The essential subtilisin-like serine protease SUB1 of Plasmodium merozoites plays a dual role in egress from and invasion into host erythrocytes. It belongs to a new generation of attractive drug targets against which specific potent inhibitors are actively searched. We characterize here the P. vivax SUB1 enzyme and show that it displays a typical auto-processing pattern and apical localization in P. vivax merozoites. To search for small PvSUB1 inhibitors, we took advantage of the similarity of SUB1 with bacterial subtilisins and generated P. vivax SUB1 three-dimensional models. The structure-based virtual screening of a large commercial chemical compounds library identified 306 virtual best hits, of which 37 were experimentally confirmed inhibitors and 5 had Ki values of <50 µM for PvSUB1. Interestingly, they belong to different chemical families. The most promising competitive inhibitor of PvSUB1 (compound 2) was equally active on PfSUB1 and displayed anti-P. falciparum and Plasmodium berghei activity in vitro and in vivo, respectively. Compound 2 inhibited the endogenous PfSUB1 as illustrated by the inhibited maturation of its natural substrate PfSERA5 and inhibited parasite egress and subsequent erythrocyte invasion. These data indicate that the strategy of in silico screening of three-dimensional models to select for virtual inhibitors combined with stringent biological validation successfully identified several inhibitors of the PvSUB1 enzyme. The most promising hit proved to be a potent cross-inhibitor of PlasmodiumSUB1, laying the groundwork for the development of a globally active small compound antimalarial.

Inhibitors of the TAM subfamily of tyrosine kinases: synthesis and biological evaluation.

The TAM subfamily of Receptor Tyrosine Kinases (RTKs) contains three human proteins of therapeutical interest, Axl, Mer, and Tyro3. Our goal was to design a type II inhibitor specific for this family, i.e. able to interact with the allosteric pocket and with the hinge region of the kinase. We report the synthesis of several series of purine analogues of BMS-777607. The structural diversity of the designed inhibitors was expected to modify the interactions formed in the binding site and consequently to modulate their selectivity profiles. The most potent inhibitor 6g exhibits Kds of 39, 42, 65 and 200 nM against Axl, Mer, Met and Tyro3 respectively. Analysis of the affinity of 6g for active and inactive forms of Abl1, an RTK protein that does not belong to the TAM subfamily, together with the binding modes of 6g predicted by docking studies, indicates that 6g displays some selectivity for the TAM family and may act as a type II inhibitor.

Identification of chemogenomic features from drug-target interaction networks using interpretable classifiers.

MOTIVATION: Drug effects are mainly caused by the interactions between drug molecules and their target proteins including primary targets and off-targets. Identification of the molecular mechanisms behind overall drug-target interactions is crucial in the drug design process. RESULTS: We develop a classifier-based approach to identify chemogenomic features (the underlying associations between drug chemical substructures and protein domains) that are involved in drug-target interaction networks. We propose a novel algorithm for extracting informative chemogenomic features by using L(1) regularized classifiers over the tensor product space of possible drug-target pairs. It is shown that the proposed method can extract a very limited number of chemogenomic features without loosing the performance of predicting drug-target interactions and the extracted features are biologically meaningful. The extracted substructure-domain association network enables us to suggest ligand chemical fragments specific for each protein domain and ligand core substructures important for a wide range of protein families. AVAILABILITY: Softwares are available at the supplemental website. CONTACT: yamanishi@bioreg.kyushu-u.ac.jp SUPPLEMENTARY INFORMATION: Datasets and all results are available at http://cbio.ensmp.fr/~yyamanishi/l1binary/ .

Relating drug-protein interaction network with drug side effects.

MOTIVATION: Identifying the emergence and underlying mechanisms of drug side effects is a challenging task in the drug development process. This underscores the importance of system-wide approaches for linking different scales of drug actions; namely drug-protein interactions (molecular scale) and side effects (phenotypic scale) toward side effect prediction for uncharacterized drugs. RESULTS: We performed a large-scale analysis to extract correlated sets of targeted proteins and side effects, based on the co-occurrence of drugs in protein-binding profiles and side effect profiles, using sparse canonical correlation analysis. The analysis of 658 drugs with the two profiles for 1368 proteins and 1339 side effects led to the extraction of 80 correlated sets. Enrichment analyses using KEGG and Gene Ontology showed that most of the correlated sets were significantly enriched with proteins that are involved in the same biological pathways, even if their molecular functions are different. This allowed for a biologically relevant interpretation regarding the relationship between drug-targeted proteins and side effects. The extracted side effects can be regarded as possible phenotypic outcomes by drugs targeting the proteins that appear in the same correlated set. The proposed method is expected to be useful for predicting potential side effects of new drug candidate compounds based on their protein-binding profiles. SUPPLEMENTARY INFORMATION: Datasets and all results are available at http://web.kuicr.kyoto-u.ac.jp/supp/smizutan/target-effect/. AVAILABILITY: Software is available at the above supplementary website. CONTACT: yamanishi@bioreg.kyushu-u.ac.jp, or goto@kuicr.kyoto-u.ac.jp.

A probabilistic model for cell population phenotyping using HCS data.

High Content Screening (HCS) platforms allow screening living cells under a wide range of experimental conditions and give access to a whole panel of cellular responses to a specific treatment. The outcome is a series of cell population images. Within these images, the heterogeneity of cellular response to the same treatment leads to a whole range of observed values for the recorded cellular features. Consequently, it is difficult to compare and interpret experiments. Moreover, the definition of phenotypic classes at a cell population level remains an open question, although this would ease experiments analyses. In the present work, we tackle these two questions. The input of the method is a series of cell population images for which segmentation and cellular phenotype classification has already been performed. We propose a probabilistic model to represent and later compare cell populations. The model is able to fully exploit the HCS-specific information: "dependence structure of population descriptors" and "within-population variability". The experiments we carried out illustrate how our model accounts for this specific information, as well as the fact that the model benefits from considering them. We underline that these features allow richer HCS data analysis than simpler methods based on single cellular feature values averaged over each well. We validate an HCS data analysis method based on control experiments. It accounts for HCS specificities that were not taken into account by previous methods but have a sound biological meaning. Biological validation of previously unknown outputs of the method constitutes a future line of work.